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Abstract 

The paper deals with homogenization of Navier-Stokes-type system 
describing electrorheologial fluid with random characteristics. Under 
non-standard growth conditions we construct the homogenized model 
and prove the convergence result. The structure of the limit equations 
is also studied. 


1 Introduction 

Rheological properties of some fluids might change essentially in the pres¬ 
ence of an electromagnetic field. For such fluids the viscous stress tensor is 
not only a nonlinear function of the deformation velocity tensor, it also de¬ 
pends on the spatial argument. A collection of interesting experimental data 
as well as a number of mathematical models of electrorheological fluids can 
be found in [TO] . 

In this work we assume that the driving electromagnetic field has a ran¬ 
dom statistically homogeneous microstructure. Then the viscous stress tensor 
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of the fluid is getting a random rapidly oscillating function of the spatial vari¬ 
ables. The corresponding system of equations takes the form (the so-called 
generalized Naviers-Stokes equations) 


( 1 ) 


du £ 
~dt ~ 
div u £ 


div^4(—, Du e )^j + div(u e ® u e ) + Vtt 
= 0, U e \dG = 0, u\ t=0 = u 0 , 


0, inGx(0,T), 


where the viscous stress tensor A(y,£) satisfies non-standard p(-)-growth con¬ 
ditions which are specified in details in the next Section. In ([TJ) u £ denotes 
the fluid velocity field and Du £ stands for its symmetrized gradient, n is 
the pressure, div(it e ( 8 ) u e ) is the nonlinear convective term, and A(x,Du £ ) is 
the viscosity stress tensor of the fluid; £ is a small positive parameter that 
characterizes the microscopic length scale. 

The goal of this work is to study the limit behaviour of u £ as £ —>■ 0. We 
assume that A(y ,£) is a symmetric matrix being a random ergodic statisti¬ 
cally homogeneous function of y G Ml In particular, the exponent p(y) that 
characterize the growth conditions of A(y,^) might be a random statistically 
homogeneous function. Under a monotonicity assumption and certain con¬ 
ditions on p , we construct the effective model and prove the homogenization 
result. We show in particular that the homogenized system is deterministic. 

Similar results in the periodic framework have been obtained in HU. 
Qualitative theory of a generalized Navier-Stokes system were developed in 
[3] and [ 12] . 

Our approach relies on a priori estimates, monotonicity arguments, gen¬ 
eralized div-curl Lemma and ergodic theorems. 


2 Problem setup 

Given a Lipschitz bounded domain G in M. d we study initial-boundary 
problem ([T]) in Qt = G x [0, T] for a fixed T > 0. 

Let (0, IF ', P) be a standard probability space with a measure preserving 
dynamical system r y , y G Ml We recall that r y is a group of measurable 
mappings r y : OgO such that 

• r yi+V 2 = T yi ° t V 2 , to = Id. 

• P (Ty(Q)) = P(Q) for any Q G IF and any y G M n . 
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• t : Q x R n K> 12 is measurable; we assume here that W l is equipped 
with the Borel cr-algebra. 

In what follows we assume that the dynamical system r. is ergodic that is any 
function which is invariant with respect to r. is equal to a constant almost 
surely (a.s.). 

We also assume that 12 is a compact metric space and that r is continuous 
with respect to this topology. 

Now we set 

where A = A(cu,£) possesses the following properties: 
hi. A : hi x IVI i —y M. where M is the space of symmetric d x d-mat rices 

d( d + 1 ) 

which is identified with R 2 . We assume that A is a Caratheodory 

function, that is A is continuous in £ for almost all uj £ hi and measur¬ 
able in u> for any £. 

h2. For all u G hi and £1 7 ^ £ 2 

( A (w £ 1 ) - A (a;, £ 2 ), 6 - £ 2 ) >0. 
h3. There exists Co > 0 such that 

(A^.O.O^ColCI^-fco)- 1 . 
h4. There exists c\ > 0 such that 

|A( W .0r V) < dlil'M + c„ p'M = /hi , 

p{uj) - 1 

where the random variable p{ui) satishes the following estimates: 

(2) 1 < a < p(u ) < (3 < 00 . 

2.1 Functional spaces 

We introduce here several functional spaces. We denote 

C„y,(G) = e C 0 “(G; H J ), div0 = 0}, 
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and H is the closure of Qj^ ol (G) in L' 2 (G ; M d ) norm. We also define X s as 
the closure of the space C' 00 ([0, T\; C^° sol (G)) in the Luxemburg norm 


\\m Lps(Q t ) 


inf |A > 0 : f |A 1 D'ip\ p ^ x ' > dxdt < 1 

1 JQt 




here Qt = G x (0,T) and p £ (x ) = p(t x / £ u>). Observe that the space X s 
depends on co. 

We say that a vector function u G X £ D L°°((0, T); H ) is a weak solution 
of problem (jT]) if 

(i) for any p G Cq^ o1 and for any t" G [0, T] the relation holds 



u(x, t')]-p(x) dx+ 



X 

A(—, Du) —u®u ■Dpdxdt 


0 ; 


n 


lim 

t -^+0 J Q 


u(x,t) ■ <p{x) dx = / u 0 (x) ■ p(x) dx 

1 Jg 

(iii) the energy inequality 

,.j" r. 

(x,t")-u(x,t")—u(x,t')-u(x,t')]dx+ / A^-,Du^j-Dudxdt< 0 


1 

2 


holds for almost all t' : t" G [0,T]. 

Notice that from the definition of a solution it follows that (u(-,t),p) is a 
continuous function of t for any p G Cq^ o1 . In other words, u(-,t) is a weakly 
continuous function of t with values in H. However, it does not imply the 
energy equality. The theory admits the strict energy inequality, which means 
the violation of energy conservation law. 

The following statement has been proved in m 

Theorem 1 Assume that 

. . d + \/3 d? + 4 d 3 d 

a > a 0 (d) = max 


d+2 'd+2 
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and 


a < p(x) < (3 < oo. 

Then generalized Navier-Stokes system m has a weak solution for any u 0 G 

H. 


Remark 1 In dimension d — 3 we have «o(3) G (1.84, 1.85). 

The condition a > a 0 ensures that the convective term u®u can be esti¬ 
mated in terms of the viscous term. More precisely, the following statement 
holds. 


Lemma 2.1 If u E X fl L°°(0, T, H), then 

\u\ 2 G L 1 (0, T, L a '(G)). 

Remark 2 In the classical case we have p = see E3i- Notice that if 
a = then 

\u\ 2 EL a \0,T,L a '(G)) = L a '(Q T ). 

In this case the convective term is completely subjected to viscous one. 

Due to Theorem [H for each c > 0 problem (jTJ) has a solution. Our goal 
is to study the limit behaviour of these solutions as e —> 0. 

The following sections deal with the homogenization procedure. This 
procedure relies on a number of auxiliary cell problems and the corresponding 
functional spaces. We introduce these spaces here. 

We denote by L p ('\Q, M d O +1 )/ 2 ) the space of functions defined on 0 with 
values in the space of d x d symmetric matrices and such that 


J dP{co) < oo. 

n 

This space is equipped with the corresponding Luxemburg norm 

= inf {a > o : J |V(co)|^"> dP(^) < l}. 

As an immediate consequence of the properties of dynamical system r and 
the Fubini theorem we have 
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Lemma 2.2 Let (p G L^(Ll, R^+ 1 )/2). Then a.s. (p{r x u) G Lf 0 ( c rxij) (M d , M^ 1 )/ 2 ) . 
Moreover, 

E / |0(r^)|^dx = |5| f \(P(u)\ p ^dP(u) 

7s 7n 

for any bounded Borel set S Clf 

We now denote by and X>j the generator of r in the z-th coordinate direction 

d 

and its domain in L 2 (f2), respectively. We also set V = f] T>i and 

i =1 

V 00 = {cfE L°°(Q) : d il ,...d ik <i>eL 2 (n) for all i,,, i k }. 


The set T>°° is dense in L p (Ll) for any p > 1. The realizations of functions 
from T>°° are a.s. smooth functions, see [5j. 

Denote Q(£l) the closure of {D^cp, (p G div^cp = 0} in 

where (D u <p)ij = \{di<pj + d)<pi), and div^d = d\(pi + ... + dd(pd- We then 
define 


g-L(n) = e L p, (-\ri;R d ( d+1)/2 ) ■. J e-vdP(u) 


0 


for all v G 



3 Homogenization 

In this section we prove a number of auxiliary statements and formulate 
the homogenization result. From item (iii) of the definition of a solution to 
problem (jT]) it follows that for each e > 0 and each w G H we have 

t 

(3) o SuP r IK|M £ (-,t)|)||^ 2(G) + f J \Du £ (x,s)\ Mx) dxds <C\\(\u 0 \)\\ 2 L 2 iG) 

0 G 

with a deterministic constant C. We recall that p e (x) = p{r x / £ uj). Consider¬ 
ing h3., h4. and (J2]) we derive from (l3|) 

Lemma 3.1 For eachco G D the sequence Du £ is bounded in L a (Qx] M fi C+ 1 )/ 2 ) _ 
and the sequence A £ = A(x/e, Du £ ) is bounded in L 13 (Q r ;M d ). 

Using the standard arguments (see [12], Section 5]), one can show that 
{ u £ (•, t)} is a family of weakly equicontinuous functions [0, T] H y L 2 (G ; M d U +1 )/ 2 ). 
Moreover, by the Aubin-Lions lemma, this family is compact in L 2 (Qt] M d ). 
This yields the following convergence result. 
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Lemma 3.2 For a subsequence, as £ ^ 0, 


u £ (-,t ) —^ u(-,t) weakly in L 2 (G;M d ) for all t G [0,T]; 
u £ (•, t) —> u(-, t) in L 2 (G;R d ) for a.a. t G [0, T]; 


Du £ Du weakly in L a (Q T ; M d(d+1) / 2 )) ; 

a(-, Diij z° weakly in l/(Q T ; M d(d+1)/2 ). 

Notice that u = u{x,t) and 2 ° = z°(x,t ) might depend on cu. 
Passing to the limit in the integral identity (i) we obtain 


(4) 



[u(x,t") 


u(x, t')\ ■ (p(x) dx + 



• Dip dxdt = 0 


for any <p G C™ sol (G) and for any t',t" G [0,T], The crucial step now is 
to determine a relation between z° and Du. To this end we consider the 
following auxiliary problem: given £ G M d ( d+1 )/ 2 find v% G Q{D) such that 


(5) 



co, v^(co) + £) • 9(co) dP(co) 


0 


for any 6 G C/(f2). 


Lemma 3.3 Under assumptions hi. h4. problem ([5jj has a unique solution 
for each £ G M d ( d+1 )/ 2 . 

Proof relies on classical result for monotone operators. Denote by A £ the 
operator mapping Q(Q) to ^(D) and defined by ^[0](o;) = A(co,f + 9(co)). 
Due to assumption h2. this operator is monotone. From h4. it follows that 
At: is bounded. Then, from hi. and h4. with the help of Lebesgue theorem 
one can derive that the function 

s — y f A (co, £ + 9\{co ) + s9 2 (co )) • 9 3 (co) dP(cu) = 0 

Jn 

is continuous in s G R for any 9\, 0 2 , 9 3 G GifTl). Also, as an immediate 
consequence of h3., we have ||0|| 1 ( v 4.^(6 ) ), 9) —y oo, as ||0|| —y oo. Then, by 
[HI Theorem 2.2.1] problem (J5]) has a unique solution. □ 

The homogenized diffusion tensor is now introduced by 

A eff (£) = / A(a;,£ + u ? M)c/PM. 

Jn 
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( 6 ) 


Consider an auxiliary variational problem 

*e0(n) J n p(v) 


The conjugate (in the sense of Young) functional takes the form 

I p'M 


no = 


\w\ 




dP(u) : w G ^ _L (r2), / wdP(uj) = 


Both functionals / and f* are convex and even. Moreover, /(£) > 0 for 
f ^ 0, and /*(0 > 0 for £ ^ 0. 

Lemma 3.4 Function f(0 satisfies the following inequality 

A“/(£), if A < 1 


/(AO < 


A /3 /(0> 


if A > 1 


Proof. Denote the function in C/(f2) that provides the minimum in 
We have 


/(AO = 


\X^ + w x ^co)\ piu) 

p(u) 


dP < 


|A£ + An; ? (a;)| p ^ 


in 


p(u) 


dP 


Jn P(v) 

This implies the desired inequality. 

Let iJ (Qt) be the associated with / Orlicz space defined as 

L f {Q T ) = \(feL 1 {Q T ,R d{d+1)/2 ) ■ [ f(<f>(x))dx< oo 

L JQt 


□ 


with the norm 


Li 


= inf \ A > 0 


/(A <f>) dx < 1 


1 Qt 


We also need the following Sobolev-Orlicz spaces: 

wP(G) = {4>€ W'/(G) : div0 = 0. f{D4>) e 0(G)}, 







and 


Xf(Q T ) = {VeL 1 (( 0,T),W 0 1 ’ 1 (G;M d )) : div x tf = 0, f(Dd) E L\Q t )}, 

= \\ m \\ L f(Q T y 

The following statement has been proved in [3} Proposition X.2.6] 

Lemma 3.5 The space C{fl sol (G) is dense in W^ (G), and the space C°°([0, T], C™ sol (G)) 
is dense in X^(Q T ). 

For star-shaped domains this result can be easily proved with the help of 
smoothing operators. For a generic Lipschitz domain the proof is more in¬ 
volved. 

The properties of homogenized diffusion tensor A eS are given in the fol¬ 
lowing statement. 

Lemma 3.6 The homogenized tensor A cS is strictly monotone and contin¬ 
uous. Moreover, the flux A(£ + vfl-)) is a weakly continuous function of £ 
with values in M d ( d+1 )/ 2 ). There exist c 0 > 0 and C\ > 0 such that 

A eS (0-Z>cof(0~Co\ 
ru !•*(«)< c t /K) + cj. 

Proof. Considering problem ([5]) and h3., we have 

dl eff (0 ■ £ = J A (u, £ + vfluj)) • £ dP = I A(u>, £ + vfluj)) ■ (£ + n ? (ca)) dP 

f2 

> Co J + vflu)\ pH dP - Co 1 > c 0 /(£) - Co 1 , 
n 

This gives the first inequality in (J7J). To justify the second one we notice 
that A(£ + v$) E and f n A(£ + v$) dP = A eff (£). Therefore, by the 

definition of /*, 

f‘(A‘ s (0) < J |A(w,C + B f (u))|f'<“>dP 

n 

< c 2 j \A(co,f + vflco)) ■ + vfluj)) dP + c 3 

= c 2 A eff (£) • £ + c 3 < c 2 ( 7 /*(A eff (0) + C7( 7 )/(0) + c 3 ; 
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here we have also used h 4 ., h 3 ., the Young inequality and Lemma 13.41 
Choosing in the last expression 7 = ( 2 c 2 ) _1 , we obtain the second estimate 

in (J 7 D- 


Strict monotonicity of A cff (£) is an immediate consequence of the strict 
monotonicity of A(ca,£) and the definition of A cS . Indeed, 

(A cff (^!)— A cff (£ 2 )H£i— € 2 ) = [ (A(u;, ^i+u^cu))—A(ca, £ 2 +^ 2 ( w )))'(£i — £2) dP 

Jn 



M)-A(w,£ 2 +w 6 (a;)))-(6+^1 (6+^2M)) > °- 


In order to prove weak continuity of A(^ + n^(-)) we first show that v^(-) is 
a weakly continuous in £ function with values in L p< d(k2, M d ). To this end we 
consider a sequence that converges to £ and notice that, due to condition 
h 3 ., we have \\v^ ||lp(-) < C. Then for a subsequence converges to some 
77 G G(Gl) weakly in L p ^(k 2 ,M d ). By monotonicity, for any ( G G(£l) it holds 


/ A(w,^+C)-(n ?J -C) dP 
Jn 



w,0+C)-A(a;,^+n ? .))-(n ? .-C)dP < 0. 


From hi. and h 4 . we deduce by the Lebesgue theorem that A (a;, £7 + £)—>■ 
A(u;,£ + <e) strongly in L p,( ^(f 2 , M d ( d+1 )/ 2 ). Passing to the limit j —» 00 in the 
last inequality yields 


/ A(w,f + CM) • (v~C(u))dP < 0 . 

Jq 

This implies with the help of Minty’s argument that 77 is a solution of problem 
( 151 ) . Since a solution of ([ 5 ]) is unique, 77 = v%. Therefore, 77 converges to v%. 

Denote by 2 a weak limit (for a subsequence) of A(-, £7+1^- (•)), as j —$■ 00. 
Since z G t/- L (Q), 


z ■ V£ (IP = 0 , / z ■ ( dP = 0 for all ( G G(Gl). 

iq. Jn 


By monotonicity, 


(A(w,fj + v^(u)) -A (w,^- + C(w))) • -CH)dP > 0 
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Passing to the limit j —> oo we get 


Z ~ A M £ + CM)) • (yz(u) - CM) dP > 0 


Using one more time Minty’s technique we conclude that z = A (a;, £T^M)- 

□ 

The homogenized problem reads 

{ c)u 

— - div(A eS (Du)) + div(u ® u) + V 7 T = 0, (x, t ) G Q T , 

div u — 0, u\dc — 0 , u\ t =o = Mo, 

We say that a vector function u G Xf(Q T ) D L°°((0,T),H) is a solution of 
problem (J8J) if 

(i) for any tp G C™ so[ (G) and for any t',t" G [0, T] it holds 

[ [u(x,t") — u(x,t')]-ip(x)dx+ f f [A eS (Du) — u®u] ■ Dp dxdt = 0; 

Jg Jt’ Jg 


lim / u(x, t) ■ <p(x) dx = / u Q (x) ■ <p(x) dx 

^+°Jg Jg 

(iii) the inequality 

- f [u(x,t'')-u(x,t")—u(x,t , )-u(x,t , )]dx+f f A eS (Du)-Du dxdt < 0 

2 Jg Jt’ Jg 

holds for almost all t', t" G [0,T]. 

We proceed with the main homogenization result of this work. 
Theorem 2 Assume that 

ad 


fj < a* = { d — a' 
Too, 


if a < d, 
if a > d 


Then almost surely, as £ —> 0, any limit point u of the family u e is a solution 
of the homogenized problem ■ 
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Remark 3 Notice that the previous theorem does not state that the limit 
function is deterministic. Although the limit problem is not random, a so¬ 
lution need not be unique. Then, the limit points of u £ might be distinct for 
different realizations. 


4 Stochastic two-scale convergence 

We first recall the definition of stochastic two-scale convergence. Let 
= v £ (x, t,uj), 0 < £ < £o} he a family of functions such that for P almost 
all uj G 12 we have v £ (-,uj) G L P (Q T ) for all £ G (0, er 0 ]. 

Definition 4.1 We say that the family v £ G L p (Qt ) weakly stochastic two- 
scale converges, as £ —> 0, to a function v = v{x,t,uj), v G L P (Q T x 12), if 

a.s. 


(9) lim sup ||u e || L p(q t) < oo, 

£—^0 

and for any <p G C™(Qt) x V 00 ^) it holds 


lim 

£— 


v £ (x, t)ip £ (x, t ) dxdt 


' Qt 


v(x, t, uj)<p(x, t, uj) dxdtdP. 


' Qt J n 


where ip £ (x,t) = (p(x,t,T x / e u). 

Notice that the two-scale limit function might also depend on the real¬ 
ization of the medium uj. Observe also that although the two-scale limit is 
defined separately for each typical realization of the medium, that is for a 
given uj, the limit function is defined on the whole 12. We do not indicate the 
dependence on uj explicitly. 


We recall some of the main properties of stochastic two-scale convergence 
(see m) that are used in the further analysis. For the reader convenience 
we provide a proof of these statements. 

Lemma 4.1 Every family of functions {v £ , £ > 0} such that (0|) holds, 
weakly two-scale converges for a subsequence to some v = v{x,t,uj), v G 
L p (Q T x 12). 
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Proof. With the help of the Birkhoff ergodic theorem we obtain that for 
any p G f> G and for almost all uj G 


lim sup 

£—>-0 


3 T 


< 


< lim sup ||u £ ||lp(q t ) ( \ \<p(x)\ q \(f)(T!Lti)\ q dx\ < 


< Ccb lim I / |v9(x)| 9 0(rEo))|' ? da: ] = C* 



\p(x)\ q \(f)(uj)\ q dP(uj)dx 1 . 


Using the diagonal procedure we can choose a subsequence e 3 —> 0 such that 
the limit lim f v £ (x)p(x)(j)(T^u})dx exists for each p and <f>. It immediately 
£ ^°Qt 

follows from the last formula that this limit defines a linear bounded func¬ 
tional on L q {Qx x f2). Therefore, there exists a function v G L p {Qt x f2) 
such that 


lim 

£—^0 


v £ (x)p(x)<f(r?LCj)dx 


Qt 



v(x, t , oj)p(x)<f{uj)dxdP. 


Qt ^ 


By the density arguments the last relation also holds for any test function 
p G Cff(Q T ) x V°°(Tl). This completes the proof. □ 

Lemma 4.2 Let a family v £ be such that a.s. 

\\v £ \\lp(q t ) < C, lim £\\V x v £ \\ L p(Qt) = 0. 

Then, for a subsequence, 

v e —^v weakly two-scale in L p (Qt), 

with v = v(x, t), v G L p (Qt). 
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Proof. Choosing a test function of the form (p(x,t)4>(r x / e oj), we get for a 
subsequence 


0 


lim 

£— 


eV x v £ (x, t)<p(x, t)(j){r x/e uj) dx 


Qt 



v £ (x , t)(p(x, t)dw u (J)(T x / e cu) dx 


v(x, t, u>)(f(x, f)div^0(a;)d£<iP. 


' Qt J £1 


Therefore, for almost all (x, t ) G Qt we have 


v(x, t, u)div u (f>(u})dP. 


in the same way as in f 131 Lemma 2.5] one can show that the set {div^d : 0 G 
T>°°} is dense in the space of L q [Ll) functions with zero average. Therefore, 
v does not depend onu. □ 


Lemma 4.3 Let a family v £ satisfy a.s the estimate 

II^IU^Qt) + || V x v £ \\lp(q t ) < C 

for all e G (0,£o]- Then, for a subsequence, 

V x v £ — 1 V X v(x, t ) + v\(x, t, c a) weakly two-scale in L p (Qt x Q), 

with v = v(x,t), v G L p ((0, T); W 1,P (G)) and V\ G L P (Q T ] L p ot (Q)), where 
Lpot(^) is the closure in L p (Ll) of the set {d^u : u G 

Proof. According to the previous Lemma a two-scale limit of v £ does 
not depend on oj. Denote by V = V(x,t,co ) the two-scale limit of V x v £ , 
and by v = v(x,t) the two-scale limit of v £ . Since the two-scale conver¬ 
gence in L p {Qt x f2) implies the weak convergence in L p (Qt), we have 
v G L p ( 0, T; W 1,P (Q)). Taking a test function ip(x, t)<f>{r x / £ uj ) with div w 0 = 0, 
we arrive at the following relation 



V (x, t, u)<p(x, t)<f>(u)dxdP 


lim 

£—^0 


V x v £ (x, t)ip(x, t)4>(r x / e u))dx 


v(x,t)'V x (p(x,t)(p(uj)dxdP = 


'V x v(x, t)(p(x, t)(j>{oj)dxdP . 


' Qt J ^ 


' Qt J ^ 
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Denoting V\{x,t,uj) = V(x,t, u>) — X7 x v(x,t) we conclude that for almost all 
(x, t ) G Qt and for any (j) G V°° such that div^ = 0 it holds 

/ V\(x, t, u)4>(u)dP = 0. 

Jn 

This implies the desired statement. □ 

Example. Periodic case. 

The periodic framework can be interpreted as a particular case of the random 
one. In this case fl = [0, l) d , J- is the Borel cr-algebra on 0, and P is the 
Lebesgue measure. The dynamical system r y is the set of shifts on the torus, 
that is for any c o G [0, l) d we set r y co = X{oj + y), where X(uj + y) G [0, l) d , 
and (c o + y) — X(u> + y) G Z d . One can observe that in the periodic case for 
any uj\ and 0J2 there exist y G such that 0J2 = r y u}\. This property plays 
a crucial role in the analysis of periodic media. 

In the periodic case Lemmas 14. 1 1 14.-II are classical and can be found in [9], 

DP- 


Considering a priori estimate (j3J) and using the arguments from [13] and 
mma, one can justify the following statement: 

Lemma 4.4 For a subsequence, 

2 

u £ ~^u(x,t) weakly two-scale in L a (Q T ), 

Du e —*■ Du(x, t ) + u\(x, t, u) weakly two-scale in L°(Qt x O), 

where ui(x,t, ■) G Q{Fl) a.a. in Qt and, 


( 10 ) 




| Du(x,t) + ui(x,t,uj)\ p ^ dxdtdPiu:) < 


00 ; 


A(-/e, Du £ ) z(x, t, uS) weakly two-scale in L 13 ( Q T x D), 


where 


\z(x,t,uj)\ p '^ 


' Qt J ^ 


dxdtdPQo) < 00, 


z(x,t,-) G ^(O) a.a. in Qt- Moreover, z 0 (x,t ) = f n z(x,t,uj) dP(co) with 
z 0 introduced in Lemma lXR 
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Proof. The two-scale convergence follows from the previous Lemmas. We 
should justify (HUD and similar estimate for z. Denote for brevity U(x, t, u>) = 
Du(x,t) + ui(x,t,u). For any 7 > 0 consider U 1 E C£°(Qt) x V°°(i}) such 
that ||U — D 7 ||i / c(Q Txr2 ) < 7. For any S E (0,1) by the convexity argument 
we have 

(11) f | (1 — 5)£/ 7 (t, x, tblj) + SDu £ (t, x) ^ T ^dxdt 

Qt 

<(1 — 5) / | U 7 ( f , x, t^uj) \ p ( T e^dxdt + 5 j \Du £ (t, x)\ p( ' T i U} ' l dxdt. 

Qt Qt 

Using the inequality \a + 5b\ p — \a\ p — Sp\a\ p ~ 2 ab = o(5) (|a| p + \b\ p ), as 5 —» 0, 
that holds uniformly in a and b , we obtain 

|(1 — 5)U 1 (t, x, teuj) + 5Du £ (t,x) | p ^ T f wl = (l — 5p(r|o;)) |U 7 (f, x, T|u;)| 

+5p(r^cu) |U 7 (f, x , T|u;) | p ^ r f 2 U 7 (f, x, r^oj)Du £ (t, x) 

+o(5)(\U' y (t,x,T2LUj)\ p ^ T e u ’^ + |Du e (t,x)| p( ' T t w) ). 

Integrating the last equality over Qt and combining the resulting relation 
with dll]) after straightforward rearrangements we obtain 

j \Du £ (t, x)\ p(j ^ dxdt > J (l — p(T|u;)) |t/ 7 (f, x, t*luj)\ p ^ t ^ dxdt 

Qt Qt 

+ f p{jiiUj)\U 1 {t,x,TiLUj)\ p<yT ^ 2 U 1 (t,x,TiLUj)Du £ {t,x)dxdt 
Qt 

+o 5 (l)(\U\t : x,r,co)\ p ^ ) + \Du £ {t 1 x)\ p ^) 1 

where 0,5(1) tends to zero as S —> 0. Due to the a priory estimates for Du e 
and by the Birkhoff theorem, the last term on the right-hand side does not 
exceed 05(1) for sufficiently small e. Applying again the Birkhof theorem we 
conclude that the first term on the right-hand side converges to the integral 



(1 - p(u)) |f/ 7 (t, x, u) \ p{uj) dxdtdP. 


Qt ^ 
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Since p{-)\W\^W can be used as a test function in the definition of two- 
scale convergence, the second term on the right-hand side converges to the 
integral 



p(cj)|Z/ 7 (f,x,cu)| p ^ 2 U' 1 {t,x,uj)U{t,x,uj) dxdtdP 


Qt ^ 


Summarizing the above relations yields 


liminf j \Du £ (t , x)| p( - T t^ dxdt > / / (1 — p(cj))|7/ 7 (t, x,a;) dxdtdP. 


Qt 


Qt ^ 



p(cj)|£/ 7 (f, x,cu)| p ^ 2 U 1 {t,x,uj)U{t,x,uj) dxdtdP + 05 ( 1 ). 


Qt ^ 


Sending first <5 —>■ 0 and choosing sufficiently small 7 > 0 we conclude that 


\UQt,x,co)\ p ^ dxdtdP < C 



Qt ^ 


with a constant C that does not depend on 7 . By the Fatou lemma this 
yields the desired statement. Moreover, we have 


liminf / \Du e {t, x)\ p>yT ^ dxdt > / / \U{t, x, oj)^^ dxdtdP. 


£—>0 


Qt 


Qt £2 


□ 


The last Lemma implies that 

(12) Du G LQQ t ), u G XQQt), Zq G Z/ {Qt)- 

Indeed, by Lemma 14.41 



f{Du) dxdt 



( min 

V w€.Q(£l) 


| Du{x,t) + w{ u;)| p(a; ) dP{co))dxdt 


< ( \Du{x,t) + ui{x,t, dP{uj))dxdt < 00 . 

J Q t 7 Jfl 7 
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Similarly, 


f*(zo)dxdt< / I \z(x, t, ui)\ v ^ dP{u)\dxdt < oo. 


' Qt 


' Qt J Q, 


It also follows from Lemma [4.41 that 


(13) 


rt 2 




zq ■ Du dxdt = 


>ti Jg 



z(Du + Mi) dxdtdP(uj). 


J gj 


Our next goal is to pass to the limit in the viscous term in (jT]). To this end 
we take the difference between the relations of items (i) and (iii) of Section 
O The resulting relation reads 


ti 



A[ —, Du £ ) ■ Du £ dxdt < 



A[ Du £ ) -u £ ®u £ 


■ V ?7 dxdt 


t 0 G to G 

l 1 \ 

[~\u £ (x,ti)\ 2 -u £ (x,ti) -r}{x)] - [~\u £ (x,t 0 )\ 2 -u £ (x,to) ■ v(%)]) dx. 


G 


for any r) G C^ ol (G). Considering the relation 
1 


(-\u £ \—u £ -r])dx = - / {\u £ ~ v\ 2 ) dx 

V2 / t=t 0 2 Jq 


ti i 

t=to 2 

and the symmetry of matrices A and u £ ®u £ , we derive 

ti ti 


t=t 0 



x 


A[—, Du £ ) • Du £ dxdt < 



x 


A[ —, Du £ ) -u £ ®u £ 


■ Dr] dxdt 


t 0 G 


t 0 G 

1 


+- / \u £ (x,to) - v(x)\ dx. 


'G 


Choosing f 0 in such a way that u £ (-,to) converges to u(-,to) in L 2 (G) and 
u(-,to) £ Wo’ a ((j), and passing to the limit e —)■ 0 yields 


1 1 


(14) 


lim 

£—>0 



A , Du £ j • Du £ dxdt < / [zq — u®u] • Dr] dxdt 


t 0 G 


to G 

+ - / \u(x,t 0 ) - r](x)\ 2 dx. 
2 JG 
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We are going to show that 77 = u(x,t 0 ) can be chosen as a test function in 
the last inequality. Let t» e a sequence of functions r] N G Qj^ ol (GO 

such that r] N —>• 77 in Wq’“(G). We substitute 77 ^ for a test function in (fT4l) 
and pass to the limit, as N —> 00 . It is clear that r] N —> u(-,t 0 ) in L 2 (G). 
Therefore, the last term on the right-hand side tends to zero. 

Regarding the convection term by Lemma 12.11 we have 

\u®u\ E L\(0,T),L a ). 



Then 

t\ t\ t\ 

[u®u)-DriN dxdt = j J \u®u) dt-DrjN dx —> J J \u®u)-Du(x,t 0 ) dxdt. 

G to G to 

By Lemma [3751 the space C^ ol (G) is dense in Wq(G). Therefore, we can 
assume that 77 ^ converges to u(-,t 0 ) in Wq(G). This yields 

ti ti 


t 0 G 



zq ■ Dr]^ dxdt 



zq ■ Du(x, to) dxdt ; 


t 0 G 

here we used the fact that 


t 0 G 


/ z 0 dt e L f *(G). 

Jt 0 

Letting t\ = to + h and combining (TT4|) with the above limit relations, we 
obtain 


lim — 

S—>0 /i 


1 1 



a(^ — , Du e ^j ■ Du £ dxdt < 


to G 


1 

h 

1 

h 
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to G 

ti 



[z 0 — u ® u] ■ Du(x, t 0 ) dxdt 
zo ■ Du dxdt — R(h) 


t 0 G 


with 


to+h 


RW = - h 



zo-(Du(x, t)—Du(x, to)) dxdt - 


t 0 G 


to+h 

'll J 

G t 0 


u®u dt-Du(x : t 0 ) dx. 
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With the help of (I13p we rearrange the last inequality as follows 
(15) 


lim — 

h 



t 0 G 




tl 

■ Du £ dxdt < - j j J z-(Du + ui) dP(uj)dxdt — R(h). 

to G fl 


Due to monotonicity of A(u;, £), for any <f> G C™(G, T>°°(Q, M d ( d + 1 )/ 2 )) we 
have 


1 1 



t 0 G 


A^-.Du^ - a(^, $(x, T x / e u)j 


■ [Du £ 


<h(a:, T x / e u:)\ dxdt > 0. 


We pass to the limit, as £ —» 0, in this relation. The term with the integrand 
A(£~ l x, Du £ ) • Du £ has been estimated in (fT5]h In other three terms we pass 
to the two-scale limit. This yields 


1 1 



t 0 g n 


A(cu, <h(x, cu))] • [Du + Mi — <3>(x, a;)] dP(u)dxdt > R{h). 


For an arbitrary Lebesgue point t 0 of functions z(-, t, •) • Du(-, t ) and z(-, t, •) • 
•) the left-hand side of the last inequality converges as / —* 0 to the 
following integral 



z(x, t 0 , cu) — A (a;, <F(x, a;))] • [Du(x, t 0 ) + u\(x, t 0 , a;) — a;)] dP(u)dx. 


g n 

It is also easy to check that both integrals in the definition of R(h) tend to 
zero, as h —> 0. Therefore, 



z(x, t 0 , cu)—A(cu, <3>(x, u))\-[Du(x, to)+ui(x, t 0 , u>)—$(x, a;)] dP{u)dx > 0 


g n 


for any test function $. By the standard Minty’s arguments 
z(x, t 0 ,u >) = A(oj, Du(x, to) + Ui(x, t 0 , a;)). 

By Lemma 14.41 we have z(x,to, •) G £/-*-(£)). Therefore, 

/ A (u, Du(x,t 0 ) + ui(x,t 0 ,uj)) ■ v(u) dP(uj) = 0 
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for any v G C/(f2), and thus Ui(x,t 0 ,oj) is a solution of problem ([5]) with 
£ = Du(x,t 0 ). We then conclude that 

z 0 (x,t 0 ) = / z(x, t 0 , dP(u) = / A(lu : Du(x,t 0 ) + ui(x,t 0 ,u)) dP(u) 
Jn Jn 

= A eS (Du(x,t 0 )). 

This completes the proof of Theorem [21 


5 Examples 

In this section we consider examples of random diffusion tensors A(x, £). 

Example 1. Voronoi-Poisson tesselation model. 

Consider a Poisson point process in with intensity 1, and construct the 
Voronoi tessellation (diagram) for this point process. It is known (see [2]) 
that a.s. the said Voronoi tessellation consists of a countable number of 
convex polytopes, we denote them H 2 ,.. .. Moreover, the polytopes can 
be enumerated in such a way that the characteristic function 1 h (y) is a 
B x ^-measurable function of y and uj for any j = 1 , 2 ,... 

Let rji, r )2 ■ ■ ■ be a family of i.i.d. random variable taking on values in 
[a,/3\ with «o(d) < a < (3 < a*. We then set 

OO 

= A (v,€) = l?l p(!,) " 1 , 

3= 1 

This diffusion matrix A = A(y ,£) satisfies all the conditions of Theorem [21 
Let (12, T, P) be the underlying probability space with an ergodic dynam¬ 
ical system r y such that p (y) = piryUj) with p(u>) = p(0). Problem ([5]) then 
takes the form: find v G such that 

j |£ + n(cu)| p( " ;) - 1 0(a;) dP 
n 


for any 9 G fy(fl). 

Taking the convolution of p with a C'^°(M d ) even function ip = p(y) such 
that p > 0 and f Rd p dy = 1, denoting the obtained function by p and letting 
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A(o;,£) = p(u) = p( 0 ), we define a diffusion matrix A(cu,£) with a.s. 

continuous in y realizations. 


Example 2. Bernoulli percolation model. 

Consider a checker board in W 1 with the cell [0, l) d . We associate to each 
cell a random variable that takes on the value 1 with probability q and the 
value 0 with probability 1 — q, and assume that these random variables are 
i.i.d. We denote these random variables by Q, j E 7L d , and the corresponding 
cells by Qj, so that Qj = [0, l) d + j. It is known (see [6]) that there is a p cr , 
0 < p cr < 1, such that for q > p cr the set \\QQj : 0 = 1} a - s - has a unique 

j 

unbounded connected component, the so-called infinite cluster. We denote 
it by C, and introduce the following two random functions: 

P (y) = a+(P- a)l c (y), a(y) = 1 + Q1 Qj (y) 

je z d 


with ao(d) < a < (3 < a*. Then 

^(2/,^) = «(2/)|e| p(y) , ^M d(d+1)/2 , 


is an admissible diffusion matrix. 
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